Signatures of Klein tunneling in disordered graphene p-n-p junctions 
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We present a method for obtaining quantum transport properties in graphene that uniquely 
combines three crucial features: microscopic treatment of charge disorder, fully quantum mechanical 
analysis of transport, and the ability to model experimentally relevant system sizes. As a pertinent 
application we study the disorder dependence of Klein tunneling dominated transport in p-n-p 
junctions. Both the resistance and the Fano factor show broad resonance peaks due to the presence 
of quasi bound states. This feature is washed out by the disorder when the mean free path becomes 
of the order of the distance between the two p-n interfaces. 



PACS numbers: 

A realistic modeling of quantum transport in 
graphene [1 has two main requirements: microscopic 
treatment of disorder and a fully quantum mechanical 
calculation of transport, taking into account the finite 
(but large with respect to the lattice constant) size of the 
samples. This is especially important close to the Dirac 
point pj, the charge neutrality point at which transport 
is dominated by evanescent modes and conventional an- 
alytical methods fail, and for interference effects. In this 
letter we report on the first such calculation in a study 
of the resistance oscillations in a disordered graphene p- 
n-p junction which incorporates both realistic Coulomb 
disorder and quantum mechanical transport. 

The physics of the p-n-p junction is governed by "Klein 
tunneling", an effect first discovered in the context of a 
relativistic particle tunneling through a potential barrier 
with a height comparable to its rest energy [3^. Since 
the occurrence of relativistic Klein tunneling is a crucial 
feature distinguishing graphene from ordinary electronic 
materials, a clear understanding of how to unambigu- 
ously observe Klein tunneling in realistic graphene sam- 
ples, where disorder is unavoidably present, is important. 
In graphene, transmission through a potential barrier has 
a pronounced angular dependence, with perfect transmis- 
sion for perpendicular incidence and a quick decrease of 
transmission probability for finite angle of incidence [4]. 
As a consequence, electrons can be confined within a sin- 
gle potential barrier, such as the p-n-p junction O |6]. 
Resonant tunneling through those confined states then 
leads to pronounced oscillations in the resistance as a 
function of system parameters El • 

The existence of this phenomenon relies heavily on the 
presence of a well defined interface between the p and n 
regions. Early experiments failed to reproduce the pre- 
dicted resistance oscillations [8 , presumably because of 
too much disorder. More recent experiments that seek 
to minimize the effects of disorder have been more suc- 
cessful [9l Uni HI] . However a theoretical understanding 



of whether or to what extent these experiments are re- 
ally observing Klein tunneling phenomena and the issue 
of the experimental conditions needed to see Klein tun- 
neling unambiguously have remained open. 

There exists strong evidence that disorder in current 
experiments is dominated by remote charged impuri- 
ties [m [13] invariably present in the graphene environ- 
ment. Due to the long range character of the Coulomb 
potential of impurity and gate charges, all scattering po- 
tentials and charge densities in the graphene sheet vary 
slowly on the scale of the lattice constant. This fact leads 
to two simplifications that make an accurate modeling of 
current experimental setups [9l HHl HI] possible. First, 
the Thomas- Fermi-Dirac (TFD) self-consistent density 
functional method is sufficient to obtain the ground state 
carrier density in the presence of disorder [HI [l5] . Sec- 
ond, intervalley coupling can be neglected for a smooth 
scattering potential, so that one can model transport us- 
ing the single valley Dirac Hamiltonian. The transfer 
matrix method of Refs. T!6^ and TT' is ideal for this pur- 
pose and can be used for experimentally relevant system 
sizes, as opposed to tight binding models which are gen- 
erally limited to small system sizes. 

Our calculation thus consists of two steps. In the first 
step the ground state carrier density n(r) and the cor- 
responding screened Coulomb potential Vgc is obtained 
using the TFD approach [14 . In the second step the re- 
sistance and Fano factor, the ratio of shot noise power 
and electrical current, are calculated using V^c as an in- 
put into the fully quantum mechanical transfer matrix 
approach [l6l [17] . We first describe the two steps of our 
method in more detail and then turn to the results. 

A schematic of the p-n-p junction under consideration 
is shown in Fig. [l^. A top gate of length Ltg is located 
a distance dtg above the graphene layer — in an air- 
bridge setup the medium in between is air with dielectric 
constant ei = 1. A back gate that controls the average 
carrier density is separated from the graphene by a Si02 
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FIG. 1: (Color online), (a) A schematic of the top gated 
graphene setup studied, (b) Profile of the screened potential 
V^c in the absence of disorder. The back gate density is fixed 
at ribg = 5 X 10"^"^cm~^ and (from bottom to top) Antg = 
-5 X lO^^cm"^ to 5 X lO^^cm"^ in steps of 2 x lO^^cm"^ 
All results in this paper were obtained for square samples of 
size W — L — 160 nm assuming Ltg — 30 nm, dtg — 10 nm, 
d — 1 nm, ei — 1 (air), and 62 = 4 (Si02). 



substrate (dielectric constant 62 =4). The effective di- 
electric constant in this setup is e = (ei+e2)/2. Together, 
the voltages on these two gates define the junction. Ex- 
amples of potential profiles are shown in Fig. ^p. 

A number of charged impurities is trapped in the sub- 
strate and just below the graphene layer. We model this 
by a random distribution C(r) of impurity charges at a 
fixed distance d with properties 

(C(r)) = 0; (C(ri)C(r2)) = n,^^5(v2 - ri); (1) 

with nimp the average 2D charge impurity density. 

In the TFD approximation, the ground state carrier 
density n(r) in the graphene layer is obtained by mini- 
mizing the energy functional 
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with respect to n(r). Here vf 
velocity for graphene, 



10^ m/s is the Fermi 



Ki(r) + Hg(r) + \j 



(3) 



with Vg = jhvpe^ Vd and V^^ the potentials induced 
by the impurity density C(r) and the top gate, respec- 
tively, and li the chemical potential. E^^\n\ is the ex- 
change correlation energy (see for details) which 
we include even though it gives only minor quantita- 
tive corrections to the transport properties of the p-n-p 
junction. The minimization is subject to the constraint 
(l/A^) n(v)£'r = nbg, where A' is the area of the sam- 
ple away from the top gate. The constraint is enforced 
self-consistent ly by varying ja. The electrostatic potential 
Vtg is expressed in terms of the top gate charge density 



Kg(r) = J 



d^r' 



^tg(rO 



(4) 



and is obtained self-consistently, for a fixed voltage dif- 
ference between the top gate and graphene, (/^tg, by re- 
quiring that in the region below the top gate Antg = 
ntg(r) — n(r) = Ctgiptg^ where Ctg is the top gate capac- 
itance. 

The potential Vsc defines a scattering problem through 
the Dirac Hamiltonian 



H = vfP • cr + Vsc{r) 



(5) 



with a = {(Jx^cTy) the Pauli matrices. The Schrodinger 
equation Hip — generates a transfer matrix A^, which 
relates the wavefunction at x = to the one at x = L, 
ipL = MipQ. In order to numerically calculate M we di- 
vide the interval (0, L) into TV equal subintervals of length 
8x = L/N and calculate the transfer matrix in each 
subinterval in the Born approximation. This gives [17| 
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where 
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Un{y) = ^ — / dxVsc{x,y). 
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We then take the limit ^ 00 in which the Born 
approximation becomes exact. From the transfer ma- 
trix Ai we calculate the matrix t of transmission ampli- 
tudes, which in turn gives us the two terminal conduc- 
tance and Fano factor as G = = {Ae^ / h)ii and 
F = iT[{l-U^)U^]/iTU^ . 

We now turn to our results. We first consider the clean 
case nimp = 0. By varying the top gate potential, the 
junction type can be varied from p-n-p to p-p-p (or n-p-n 
to n-n-n depending on the back gate voltage), as shown 
m Fig.[l]D. The depend ence of the resistance R and Fano 
factor F on the two gate voltages is shown in Fig. [2^-b. 
In the p-p-p {n-n-n) region of the upper left (lower right) 
corner the resistance is small compared to the resistance 
of the p-n-p {n-p-n) region of the upper right (lower left) 
corner. Pronounced oscillations in the resistance are seen 
in the p-n-p and n-p-n regions. The behavior of the Fano 
factor mostly follows that of the resistance. 

The resistance oscillations can be understood as aris- 
ing from resonances through quasi bound states inside 
the barrier created by the two p-n interfaces O [H [7]. 
The larger the transverse momentum Qy, the larger the 
incident angle for scattering off the barrier and the tun- 
neling amplitude decreases. The broad oscillations arise 
from resonant tunneling of the few modes with smallest 
Qy. Larger Qy give rise to very narrow resonance seen in 
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FIG. 2: (Color online). Gate voltage dependence of (a) re- 
sistance and (b) Fano factor of a clean device. A single trace 
of (c) resistance and (d) Fano factor at a fixed back gate den- 
sity ribg = 5 X lO^^cm"^. The smooth curves are obtained by 
averaging over boundary conditions. The insets show a close 
up of one of the narrow resonances obtained for transverse 
periodic boundary condition. 



Fig.|2]3-d, where we show a cross-section of Fig.[2^-b at a 
fixed value of the back gate voltage in a higher resolution. 
Since the positions and widths of the narrow resonances 
are sensitive to the transverse boundary conditions we 
also plot the smooth curves obtained by averaging over 
twisted boundary conditions. The weak oscillations in 
the p-p-p and n-n-n regions and the narrow resonances 
are absent after the averaging. Observation of the nar- 
row resonances will thus require very well defined edges. 
The double peak structure of the Fano factor resonances 
reflects the increase of the transmission probability from 
zero to one and back to zero and the fact that the Fano 
factor of a perfectly transmitted mode is zero. 

In contrast to the narrow resonances the resistance 
minima for the broader oscillations do not depend on 
boundary conditions and can be estimated by a semiclas- 
sical argument [5 . One can think of the p-n-p junction as 
a Fabry-Perot etalon in which waves scattered off the two 
p-n interfaces interfere destructively if the WKB phase 

OwKB = - Vsc{x', y)dx' (mod tt), (8) 

J Xi 

with the turning points Xi defined by the condition 
Vsci^ i^y) = ^ = 1^2. The positions of the resistance 
minima that follow from this argument are indicated by 
arrows in Fig. [2] 

We now consider the effect of disorder. A single re- 
alization of the scattering potential V^c in a disordered 
p-n-p junction is shown in Fig. [3^. For realistic impurity 
concentrations, the correlation length of V^c is ~ 10 nm 




FIG. 3: (Color online), (a) Density n(r) for a disordered 
junction with ribg = 5 x 10"^"^cm~^, Antg = 2 x 10"^^cm~^, 
and riimp = 5 x lO^^cm"^. (b) resistance and (c) Fano factor 
in the presence of disorder (single realization) for a fixed back 
gate density ribg = 5 x lO^^cm"^ and several values of the 
impurity density (from bottom to top riimp = 0,1,2.5,5,10, 
and 15 x lO^^cm"^). The disordered curves have been shifted 
for clarity, with dashed lines showing their zero, (d) The 
top gate density dependence of the resistance at a narrow 
resonance for different values of the impurity density. 



[m [19] . While one can still clearly make out the different 
p and n regions in the presence of disorder, the boundary 
between the two is no longer as sharp, leading to a weaker 
confinement of particles inside the barrier. The effects of 
Klein tunneling are thus expected to be suppressed. 

Figures ^jp-c give the transport properties as a func- 
tion of top gate voltage for a single disorder realization 
and for different impurity densities. Mesoscopic fluctua- 
tions due to the disorder are superimposed on the oscilla- 
tions from the Klein tunneling. As the impurity strength 
increases these fluctuations become stronger, eventually 
dominating the signal. The Fano factor is more sensitive 
to disorder than the resistance. The effect of disorder on 
a narrow resonance is shown in Fig. [3]i. 

In Fig. [4] we show the disorder averaged resistance and 
Fano factor, where the average was taken over 10^ dis- 
order realizations. We find that the clean limit broad 
quantum oscillations survive for impurity densities up to 
^imp ~ 10^^ cm~^ while the narrow resonances disappear 
at nimp ~ 10^^ cm~^. We do emphasize that for weak 
enough disorder the narrow resonances will be present, 
and it is therefore conceivable that they may actually be 
observed at low enough temperatures, when the phase 
coherence length is much larger than the system size. 

As in the case of the single disorder realization, the 
dependence of the Fano factor on the disorder strength 
is similar to that of the resistance. The dependence of the 
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tunneling, to the disordered diffusive transport is found 
to take place as the mean free path becomes of the order 
of the distance between the two p-n interfaces consistent 
with recent experiments . 
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have been performed on the University of Maryland High 
Performance Computing Cluster (HPCC). This work is 
supported by US-ONR and NSF-NRI, by the NSF under 
grant no. DMR 0705476, and by the Humboldt Founda- 
tion. JHB thanks the Dahlem Center at FU Berlin for 
hospitality. 



FIG. 4: (Color online). Disorder averaged resistance (a) and 
Fano factor (b) as a function of top gate voltage for different 
values of impurity densities (parameters and color coding the 
same as in Fig. [s] b. (c) Disorder averaged resistance at a 
narrow resonance for different values of the impurity density, 
(d) Top-gate density dependence of the disorder averaged re- 
sistance for different values of back gate density. 



resistance on the back gate voltage for a fixed impurity 
density is shown in Fig. [4]i. The main effect of the back 
gate is an overall reduction of the resistance and a shift of 
the top gate density at which the central region changes 
polarity. 

One expects impurity scattering effects to start domi- 
nating the ballistic oscillations when the mean free path 
/ becomes comparable to the sample width W (for the 
narrow resonances) or the length of the central n region 
(for the broad oscillations). The mean free path / can 
be estimated through its relation to the conductivity a, 
/ = {h/2e^){cF /kp)^ where /ci? is the Fermi wavevector. 
Away from the Dirac point, a = {2e^ /h){n) /ni^^f{rs^d) 
[20] , where /(rg, d) is a known function of Vg and d only 
[f{rs = 0.8, d= Inm) = 0.1] [22^. Setting (n) = nbg and 
using the impurity densities quoted above, we find that 
the narrow resonances vanish for / ^ 400 nm, which is 
the same order of magnitude as expected. The impurity 
density at which the broad oscillations disappear corre- 
sponds to a mean free path / ^ 40 nm. This is roughly 
the same as the length between the two p-n interfaces, 
see Fig. [l] and consistent with the experimental results 
of Kefs. H and [lOl 

In summary, we have presented a powerful theoretical 
method that is generally applicable to transport prob- 
lems in realistic graphene samples where the disorder 
is dominated by charge impurities and transport prop- 
erties need to be obtained fully quantum mechanically. 
We have applied this method to understand the effects 
of disorder on transport through p-n-p junctions. The 
crossover from the ballistic transport governed by Klein 
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